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We compute cross sections for inclusive scattering of high energy electrons on ^He, based on the 
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density matrices are obtained from non-relativistic *He wave functions using realistic models for the 
nucleon-nucleon and three-nucleon interaction. Predictions for E=3.6 GeV agree well with the NE3 
SLAG- Virginia data. 
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> ■ I. INTRODUCTION 

o ■ 

, Total nuclear structure functions (SF) contain at least two components, namely one describing nuclei as composed 
\^ ' of point-particles and a second one accounting for the internal structure of the nucleons. Since the latter are taken 
, from experiment, a computation of total nuclear SF amounts to a determination of the same for a nucleus which is 
composed of point-particles. Such a calculation is usually performed within one of the following two approaches. In 
' the first, one perturbatively expands the SF in the residual interaction between a nucleon struck by the virtual photon 
, and the remaining spectator nucleus, thus generating the Impulse Series (IS) for the SF. In the widely used lowest 
' order Impulse Approximation this residual interaction is first neglected. One then either computes higher order Final 
[ State Interaction (FSI) terms (see for instance Refs. ||l|,||), or models them [||-|6|. 

An alternative approach is based on a relativistic generalization of the Gersch-Rodriguez-Smith (GRS) expansion 
of SF in inverse powers of the 3-momentum transfer |<7| [0-^. Both theories have been applied to cross sections for 
inclusive scattering of high-energy leptons from various nuclear targets [PH6 HlO|JTl|l . 

■ When applied to high energy inclusive scattering one usually limits a GRS calculation to the two lowest order 

■ terms. Their determination requires knowledge of one- and two-particle density matrices, which are not diagonal in 
the coordinate of the struck nucleon '1', and of the spectral function. The non-diagonal one-body density matrix 
is related to the single-nucleon momentum distribution n{p) and is usually extracted from alternative experimental 

^ sources, or is computed from theoretical models. There generally is no direct information on the half-diagonal, two- 
(— I particle density matrices for finite systems and one relies on parametrizations [^,p^. In those, nuclear recoil is usually 
neglected, thereby limiting applications to targets with A > 12. 

In the following we exploit accurately computed non-relativistic (NR) wave functions for light nuclei, using a number 
of modern realistic nucleon-nucleon {NN) and three-nucleon (3A^) interactions. Those wave functions are Galilean 
^ ■ invariant and enable a realistic GRS calculation of inclusive scattering on those nuclei. As a first application we 
5^ , choose "^He and for that target we shall report below predictions and a comparison with data. 

The other ingredient, namely the "^He spectral function V{p, £), is rather difficult to compute and only a few direct 
calculations are reported Below we shall adopt the reasonable alternative that has been developed by Ciofi 

degh Atti and Simula 

At this point we mention that for years the IS and GRS approaches have been considered as being distinct and even 
incompatible. Only recently has their equivalence been demonstrated, provided both series are expanded to the same 
order in the same parameter Following the derived prescription to link the two approaches, one can perform 

an interesting numerical comparison. 

The GRS and IS theories are not the only tools which have been used to compute nuclear SF. We mention in 
particular the ingenious method of Efros and co-workers, which has been applied to "'He Regrettably, it 

appears not feasible to extend that method to high energies. 

The present note is organized as follows. In section O we recall the GRS approach, emphasizing the two main 
ingredients of our calculations, namely the SF of a target composed of point-particles and the SF of the free nucleons. 
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We also discuss there the computation of the above density matrices. In section [II we compare predictions for cross 



sections with the Virginia-SLAC data M. In the last section we present our conclusions. 
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II. TOTAL NUCLEAR STRUCTURE FUNCTIONS. 



The cross section per nucleon for inclusive scattering of high-energy electrons from a nucleus with A nucleons reads 



d^aeAiE;e,iy)/A ^ 2_ 
dn dv M 
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(2.1) 



where M is the nucleon mass, aM the Mott cross section, E the beam energy, 9 the laboratory scattering angle and 
ly the energy loss imparted onto the target. The above nuclear structure functions per nucleon F^{x,Q^) contain 
the essence of unpolarized electron scattering from randomly oriented targets. Those SF depend on the squared 
4-momentum transfer —q^ = = — i/^ and on the Bjorken variable x = /2Mv with range Q < x < A. For 
given beam energy i?, [9, v) and (x, Q^) are sets of alternative kinematic variables. 

Total nuclear structure functions per nucleon may, in a semi-heuristic fashion, be expressed as follows [p^-p^ 

•'^ 1=1 ^ ^ 

where fP'^'^ is the SF of a nucleus composed of point-particles and Plf^'^ is the averaged free nucleon SF. For a 
nucleus A{Z,N) 

i^f > {x, Q') = jFJ:{x, Q') + ^F-{x, Q2) , = 1, 2 . (2.3) 

The coeffici ent functions Cki account for the mixing and modification of the free nucleon structure functions in the 



expression ( p. 2]) [pO[ . We retained in this paper only the dominant coefficient 



C22{^,Z) « [1 - p + p/z? - ip(l - p)(l - l/zf , 

p= [l + Q'^/AM^x^]-'^ . (2.4) 

Eq. (p.4|) is a better approximation for C22 than previously used . 

Eq. ( |2.2| ) is valid for x > 0.15 — 0.20, below which pionic and anti-screening effects become of importance |^l| and 
above some critical Q^, which presumably can be estimated using QCD. A previous comparison of predictions and 
data for medium- A targets produced an empirical estimate ^ Qc ~ 2.0 — 2.5 GeV^ |11[. 



Each of the SF F^'" in (2.2) has both nucleon-elastic (NE) and nucleon-inclastic (NI) parts, thus F^ = F^. 



F^'^^, with N = p,n. The total nuclear structure functions, Eq. ( |2.2D , and the total cross section per nucleon may 
therefore be expressed as a sum over contributions coming from the NE and NI parts of nucleon SF. In particular, 
the NE part p^'^^ is the well-known combinations of static electro-magnetic form factors and contributes primarily 
around the region of the quasi-elastic peak (QEP), x ~ 1. For the inelastic parts Fj^'^^ we have taken 

F^''^^(x,g2)«F,^-^i(x,Q2)-Ff^^(x,Q2) ; fc = l,2, (2.5) 

where F^'^^{x, Q^) are the deuteron SF's per nucleon. For F^'^^{x, Q^) and F^'^^{x, Q^) we employ values interpolated 
between the data of Ref. ||2^, whereas for ^^'^^(x, Q^) and F^'^'(x, Q^) we use the parametrizations of Ref. Q. 

A. The GRS series 

We now focus on f^^'^ in Eq. ( |2.2| ), the SF for a nucleus of point-particles, which has to be computed. Following 
Ref. @ one writes 

fP^'A{x,Q^)^(^) 0(|gUG), \q\ = Q^l + {Ql2Mxf , (2.6) 

\ ox J Q2 fixed 

where (j){\q\,y) is the reduced response in terms of a relativistic scaling variable H 
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(2.7) 



and (A) some a vera ge nucleon separation energy. We shall retain the above 1 /A' correction in the scaling variable yc 
which, as Eq. ( ^.7|) shows, is simply related to the Bjorken variable x. In the GRS approach the reduced response 
may, for smooth NN interactions, be expanded in a series of inverse powers of |q| Explicitly 



The lowest order term is given by 



|gUG°) = 



dp p 
4^ 



d£P{p,£) + e{y^") 



Vg 



dp p 
4^ 



d£V{p,S) , 



(2.8) 



(2.9) 



with V{p,£) the standard single-hole spectral function. The energy argument \s £ — E — with E the removal 
energy and Aq the (p, n averaged) minimal separation energy (for ^He Aq s 
variable give in Eq. (EJ) with (A) = Aq. The integration limits in (^.9[) are 



20.2 MeV). Above, is the scaling 



Em {ycP, \q\ 



{yG±p)\q\ 



In actual calculations the spectral function has been written as in Ref. |2j 

r{p,£)=no{p)S{£)+ri{p,£) , 



(2.10) 



(2.11) 



where no{p) is the partial momentum distribution due to intermediate states of one nucleon and the A — 1 spectator 
system in its ground state. Contributions from continuum states of that system are summed in Vii^p, £). As stated in 
the Introduction, that part of the spectral function for ^He has been taken to be Vi{p,£) = N'{p)V^°^°^{p,£), where 
V™°'^°^{p, £) has been provided by us by C. Ciofi degli Atti 0. The normalization factor A/'(p) is fixed by 



d£ Viip,£) = n{p) - noip) , 



(2.12) 



where the quantities n{p), the total mo ment um distribution, and no{p) have been calculated using the NR wave 
functions as will be explained in Subsect. II B . 

We have also tested the following simple 2-state approximation for the spectral function [p3jlq| 



r{p, £) « m{p)6{£) + [n{p) ~ no{pm£ - (A) + Aq) , (A) « 50 MeV . 
Substitution into (|2.9|) produces a |<7| -independent lowest order contribution. 



Since in the relevant p region ^ n, an accurate approximation of Eq. (2.14) reads 

<t>^o\yG)^Y^ / , dp p nip) . 
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(2.13) 



(2.14) 



(2.15) 



Terms with n > 1 in Eq. (2.8) describe FSI corrections to the asymptotic limit as a series in l/|q|. It is easier to 
give those in terms of their Fourier transform 0„(s), namely 



Mvg) 



ds 
2^ 



0n(s) 



(2.16) 
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Each (finis) = 4>n{s; [V]) is a functional of the bare interaction V, for instance 

dri J dr2 P2iri,r2;ri ~ sq,r2)[iXqib, z; s)] , 



M 



I daV{b,z-a) 
Jo 



(6, ^; s) = ^^sVib, z-s)^ 



(2.17a) 
(2.17b) 
(2.17c) 

(2.17d) 



where b (z) is the component of the vector r = r i — r2 perpendicular (parallel) to the q direction, and p2 the semi- 
diagonal two-particle density matrix. Eqs. ( 2.17 ) define two parts of the off-shell eikonal phase Xq which are related, 
and thus 



XqK 



d 



{b,z;s) = [1- s— ]xl^\b,z,s 



(2.18) 



One frequently deals with interactions V which have a strong short-range repulsion (or produce for other reasons a 
diffractive elastic amplitude) and it is then of advantage to perform a summation over a ladder of bare interactions V. 
The replacement V — > Veff = tq, produces a well-behaved, g-dependent, off-shell t-matrix as an effective interaction, 
which in coordinate space is proportional to the off-shell profile function f ||l^. For the part f generated by 
one has 



tx'i'' (b, z;s)-. f « (b, z;s) = e^^x'^^ (b, z;s)]-l 
«0(z)0(.-z)rW(6). 



(2.19a) 
(2.19b) 



The approximation Eq. (2.19b) has been tested in Ref. ||25|. Its application permits the exploitation of a standard 

parametrization of the on-shell profile rq^\b) in terms of elastic NN scattering data, as are aq°*,Tq,Q'^\ which are 
respectively the total cross section, the ratio of the real to imaginary part of the forward elastic amplitude and the 
width of the diffractive amplitude. Explicitly, 



l(0)l2 



An 



(2.20) 



T here i s no simple way to generalize ( |2.18| ) to the total off-shell phase x- Yet as in [|l^ we shall assume that 
Eq. (2.18) is also approximately valid for the total off-shell profile function 



-iXqib, 



z: s 



tqib,z;s) 



1 - s 

1 
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fW(b,z;s) , 

e{z)eis-z)TW(b) 



(2.21a) 
(2.21b) 



After substitution of the above expression in Eq. (2.17), the leading FSI contribution to 0(q, s) turns into the following 
q-dependent result ||l0|| 



where 



M ~ , , M ~ , 



Gii\qls)^G['\\q\,s)+GY'i\q\,s) , 



(2), 



dri I dr2 P2{ri,r2;ri - sq,r2)e{z)0{s ^ z)Tip{b) 



Gf \\q\,s)^ J dri J dr2 p2(ri, ra; ri - sq, r2)s0(z)5(s - z)rW(6) . 



(2.22) 

(2.23a) 
(2.23b) 



Previous analyses dealt with targets with A > 12. For those there do not exist computations from first principles 
for single-nucleon momentum distributions uqIp), n{p) and density matrices p2, as required in ( 2.11| ), (2.12), (2.23). 
Moreover, suggested parametrizations |0,|2l do not account for nucleon-recoil, which is only justified for A > 12. 
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One of the predictions from Eq. ( p. 2D is a weak-^ dependence of the SF for point-nucleon nuclei and of the averaged 
nucleon SF |lC|]. This entails predicted inclusive cross sections per nucleon to be practically independent of A. 
Supporting evidence comes from experimental ratios of cross sections per nucleon for different targets at identical 
kinematical conditions . Definitely larger deviations from smooth A-dependence are expected, if one of the 

targ ets is a light nucleus with A <6. This is evident from Table || where we entered some C/Fe ratios from JLab data 
13 and for Hc/C from the older NE3 data Ref. 

For the above reasons we did not include in the past a GRS analysis of inclusive scattering on the lightest targets. 
In the following we exploit the possibility to compute a precise NR nuclear ground state wave function $0 of light 
nuclei f or giv en nuclear interaction V. Those enable a calculation of n{p), n(){p) and /92, which enter the components 
(PH), dill) and (|^) of the nuclear SF. 



B. The density matrices 

Various methods permit nowadays an accurate calculation of the ^He ground state wave function |^^. We exploit 
here the Correlated Hyperspherical Harmonic function (CHH) technique which has been developed by the Pisa group. 
The spatial configuration of the system is described in terms of a given choice of the Jacobi vectors ^1, ^3- the 
hyperspherical framework we use as new variables the hyperradius p, defined by 



(2.24) 



and the set Q — {^1,^2, £.3,^2,^3}- The latter includes the polar angles £_i = {Oi,(t)i) of each Jacobi vector and 
additional hyperspherical angles (p2, f3- We then write for the ground state wave function $4 



$4 



Ntot 

E 

n=l 



Un{p) 



A\ J'>,(ri2,ri3,ri4,r23,r24,r34) 3^„(il) 



(2.25) 



where A is an anti-symmetrizer and yn{^) are the four-body Hyperspherical Harmonic (HH) functions | |29| ]. Choosing 
J-ji = 1, Eq. ( 2.25 ) generates an uncorrelated HH expansion for the ^He ground state wave function. For it, the rate 
of convergence is extremely slow when the NN interaction is strongly repulsive at small distances. One accounts for 
the latter property by multiplying every HH function in the expansion with a suitably chosen correlation factor J-'n, 
ultimately leading to the CHH expansion. The latter much improves the description of the target wave function for 
small inter-nucleon distances, and a much smaller number of basis functions is required toget convergence. 
In the case of ^He, the correlation factors have been chosen to be of the Jastrow form pQ] 



J^n = fn{ri2)9n{ri3)9n{ri4.)9n{r23)9n{r24)9nir34:) 



(2.26) 



i.e. products of one-dimensional functions /„ and 9n which are solutions of a NN Schrodinger-like equation (for 
details, see Ref. | 
Using the Ray^ 



eigh-Ritz variational principle for varying functions m„ in (2.25) 



< Su<^>4\H - E\^4 >=0 



(2.27) 



one is led to a set of hyperradial equations for the functions u„ in the variable p which, after discretization, is converted 
into a generalized eigenvalue problem and are solved by standard numerical techniques . One thus determines the 
hyperradial functions Un{p) in Eq. (2.25) and the bound state energy E. We shall present below results based on 

1) the Argonne V18 nucleon-nucleon (NN) potential (the AV18 model), 

2) the Argonne V18 NN potential supplemented by the Urbana IX three-nucleon potential (the AV18UR 
model) , 

3) the Argonne V14 NN potential Q plus the Urbana VIII three-nucleon potential ^ (the AV14UR model). 

The two models, which contain a three-nucleon interaction provide a *He binding energy rather close to the experi- 
mental one, whereas the AV18 under-binds by about 4 MeV. Th e pres ent status of the ^He binding energy calculations 
with the CHH method is summarized in Table ||, where in Eq. ( ^.25 ) up to Ntot ~ 200 functions have been used (the 
explicit CHH states included in the expansion are discussed in Ref. |]30[|). 

Calculated binding energies for the AV18 or AV18UR Hamiltonians are within 1 % of the "exact" Green's function 
Monte-Carlo (GFMC) results j^] for corresponding interactions. Somewhat less satisfactory agreement between the 
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CHH and GFMC results for the AV14UR model, since this interaction is more repulsive at short distances than the 
other two. For all we checked that our final results for the deep inelastic scattering cross sections depend only slightly 
on the value of Ntot, once Ntot ^ 200. 

The thus constructed ground state wave function, readily gives the corresponding 4-body density matrices, in 
particular the one non-diagonal in '1' 



/04(ri, r2, ra, r^] r[,r2,r3, ri) ^ $4(^1, ^2, ra, r4)$4(r'^, r2, rs, r4) 



(2.28) 



Successive integrations over the diagonal coordinates 3,4, and eventually over coordinate 2, then furnish p2 and pi 
non-diagonal in '1'. The total momentum distribution is the Fourier transform of pi. 

The partial momentum distribution no{p) is obtained from the overlap of $4 and the (^H, ^He averaged) three- 
nucleon ground state wave function $3, namely 



"-o(p) = Hp)? 



a{p) = ( $3(1, 2, 3)xAm\jo{pRi23A)M^, 2, 3, 4) ) 



(2.29) 



where Xi iVi) is the spin (isospin) state of particle 4, jo the zero-order Besscl function and i?i23,4 the distance of 
particle 4 with respect to the center of mass of the other three. The ground state wave function $3 of the three- 
nucleon system has been obtained with the s ame Hamiltonian model used to generate $4, and again by application 
of the CHH technique In equation ( 2.2E ), $3 and X4^4 are coupled to give a state with vanishing total angular 
moment um a nd isospin. 

Eqs. ( 2.23| ) for G invite to express p2 in terms of the variables r — ri ~ r2 and R — (ri -I- r2)/2 and then to 
perform the R integration 



B{b,z;s)^ j dR P2{ri,r2;ri - sq,r2) 



(2.30) 



The result is substituted into (2.25) and subsequently integrated over b and z. 



C. An effective IS series 



In Ref. m one may find a detailed account of the lowest order PWIA calculations as part of the IS series. To it one 
should perturbatively add FSI due to the interaction of the knockcd-out nucleon and the spectator core. 

In the Introduction we recalled a recent proof that the GRS and IS, both evaluated to order C'(l/|<7|^), produce 
the same result. The lowest order is given by 



'/'o(|g|,2/o) = 

where in this case 



m\+yo 



\yo\ 



dp p 
4^ 



d£Vip,£)+9iy„ 



dp p 
4^ 



d£ V{p,£) 



Em = v/M^ + (yo + 
Above, yo is the IS scaling variable defined by |^,D 



(2.31) 



(2.32) 



Vo 
So 



-\q\ + ^{iy-Ao)iiy~Ao + 2M) , 

1 
2A' 



1 + v/M 



(2.33) 



.i + |g|/yo°° 

The IS first order is then obtained f rom the prescriptio n whi ch puts G^^^ = in ( p. 22 ) and simultaneously replaces 



the GRS scaling variable Eq. (^.7D, by j/o in Eq. ( 2.16 ). Finally, f^'^'^{x, Q^) is now calculated using again 
2/0 in place of yQ in Eq. (p 6[). Those changes have been shown to generate the IS series to order l/|q| Actual 
applications will be given in the following section. 
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III. RESULTS 



In this section we present numerical results for inclusive cross sections on '*He based on Eqs. (2.1), (2.2). Unless 
stated differently, t hose results use density matrices based on the AV18UR interaction and the full spectral function 
as discussed in Sec. II A. In the GRS approach, the cross sections have been computed by approximating the reduced 



structure function as follows 



order , 
0+1 order 



(3.1) 
(3.2) 



where (^o is given by Eq. (2.9) and Gi is the inverse Fourier tran sform of the functions given in Eq. ( |2^ ). The 



corresponding IS expressions are obtained as discussed in Sect, p C| . 

In Fig. 1 we display the SLAG- Virginia NEB cross section data jl^] for i?=3.595 GeV, scattering angles 9 =16°, 
20°, 25°, 30° and varying energy loss and also our numerical results obtained with the above reduced response 0oi- 
Table ID shows the ranges of the kinematical variables x,v,Q^. For all measured scattering angle, and in addition 
for 9 = 45°, we entered there values of the energy loss v and of close to the appropriate lower (x — 2) and upper 
{x — 0.1) ends of the theoretical curves shown in Fig. 1. We also show the v, values corresponding to the position 
of the QEP at X = 1 . 

Let us first focus on the NE part of the cross section (thin dashed lines), which contribute primarily around the 
QEP. As can be seen from the figure, the NE parts well describe the QEP, in particular at low . The data for 
9 = 16° show a clear maximum and an adjacent minimum which get fuzzier and ultimately disappear for increasing 
9 or . The same maximum occurs in the inclusive cross sections on deuterons [^^, but is absent for targets with 
A>12 0,13. 

This structure is the result of the competition between the NE and the NI parts of the cross section. For decreasing 
the NE part around the QEP x = 1 grows relative to the, usually dominant, NI part. For the smallest in the 
data, the NE part beyond the QEP stands out until for a; < 1 the NI one overtakes. It is the maximum value of 
jPN,A(^^^ Q2-J 2; ~ 1) which sets the magnitude of the QEP. For a given Q^, that peak value decreases with A: it 
is largest for the deuteron and "'He and then it is almost independent on Aiox A> 12, reflecting the smearing of the 
momentum distribution due to the Fermi motion. For example, f^^'^{x — 1, Q^) for the deuteron and ''He is « 5.5 
and 2.2 times, respectively, larger than for nuclei with A > 12. It causes the QEP for equal kinematic conditions to 
be most prominent for the lightest nuc lei. 

We already mentioned that Eq. (2.2) (for the NI part) is estimated to be valid for x > 0.2 and above some critical 
« 2.0 — 2.5 GeV^ [|l^,|ll[. As Table III shows, that approximate critical value is actually never reached for any 
9 = 16°, 20° data point, which renders those data not really suitable for a test of the theory. For the same reason we 
excluded from our analysis NEB data at lower energies and the same is the case for old, near-elastic, high-Q^ data 
on ^He |3^. For both angles, the convolution formula predicts too large NI parts around the QEP, resulting in an 
over prediction of the data in that region. For 9 = 30°, on the other hand, > and, in fact, comparison of data 
and predictions shows that there is good agreement for all but the smallest energy loss values. Since cross sections 
there have fallen by orders of magnitude, one expects s ensi tivity to small dynamical details. For example, without 
the inclusion in Eq. (2^) of the mixing factor C22, Eq. (2^), the agreement would be of definitely lower quality (see 
also below). 

A more stringent test for the theory would be provided by data at higher beam energies with, in general, higher 
Q^. Unfortunately, the recent 4 GeV experiments at JLab for various targets did not contain ^He but a recent 
JLab proposal includes that target in a 6 GeV run with scattering angles 9 — 15°, 23°, 30°, 45° and 60° ^Sj . The 
kinematical region explored by that experiment covers 0.2 < x < 1.0 and 1 < < 8.0 GeV^ (Table 1rv|) and 
predictions for the four largest scattering angles can be found in Fig. 2. Incidentally, we checked that for £' = 6 GeV 
the effect of C22 is practically negligible due to Q^, which grows with beam energy E. 

Next we discuss the effec t of t he diff erent approximations for 0o presented in Sec. ||. The cross sections calculated 
using the approximations ( 2.14 ) and ( 2.15|) fo r (j)o almost overlap. Moreover, they are rather similar to the ones 
computed using the spectral function, Eq. ( ^.9[ ), except in the low z/ region. There, the use of the "full" model slightly 
reduces the cross sections. 

This result deserves some comment. The differences between the two expressions, Eqs. (2.9) and ( ^.15 ) are generally 
sizeable in particular in the low v region |^ where yQ° is negative and large in absolute value. However, in the 
kinematical region of the NEB experiment, the values of Ei\i entering Eq. (2.9) are found to be rather large and then 



d£V{p,£) 



d£V{p,£) = n{p) 



(3.3) 
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As a result, the (f>o calculated using Eqs. ( p.9[ ) or (2.15) nearly coincide. For example, for E = 3.595 GeV, 9 = 30° 
and 1/ = 0.7 GeV, ~ 3 and Em is large for all the values of p. For the same reason we expect that the predicted 

cross sections do not much depend on the parametrization chosen for the spectral function, in particular not in the 
low h' region. In that region they are rather sensitive to the tail of the momentum distributions, which in turn is 
related to the correlations in the nuclear wave functions . 

In Fig. 3a we display the separate contributions of the cross sections at E — 3.595 GeV, 6 = 30° and varying v, 

as due to the NE and NI components of the nuclcon SF fI^^K The thin (heavy) dashes are NE part in the 4>o (0oi) 
approximation for the reduced response. Those have their maximum at x « 1, and are only in the wings marginally 
affected by the 1st order FSI terms. The thin and heavy solid lines show the corresponding NI parts, which by nature 
dominate the region x 1 for relatively high v. FSI affect only the low v region and cause a rather small increase 
in cross sections. In those i/-regions NE and NI parts are of the same order. Fig. 3b is as Fig. 3a for 9 — 45° (this 
angle was chosen since the values are larger and similar to those for E — 6 GeV, 9 = 15°). The results are similar, 
except that FSI now appear to decrease the 0th NI contribution at low ly. 

From Fig. 3a we observe that the slight over-prediction in the low i/ region of the theoretical results a.t 9 — 30° is 
mainly due to the NI part and as discussed above, is only marginally affected by the inclusion of FSI. The reason of 
the large NI contribution can be simply understood by looking at the convolution (E^) between the nucleon SF and 

jPN,A 

Typical behavior of the functions f^^'^{z, Q^), F^^'^'^^ {x/ z, Q^) and the mixing factor Ci^iQ'^ I x^ ^ z), (the latter 
two for x = 0.1,2.0) are given in Fig. 4 [p^'^'^^ behaves similarly). For x < 1, the permitted range of values 
z^x covers the z « 1 region where f^^'^{z^ Q^) is large and allows virtually the entire support x/ z of to 
contribute. In contrast, for a; > 1 (see Fig. 4b), only the tail of f^^^^ contributes and i^^^^'^^(i, Q^) for x/A < t < 1 
is usually small. Moreover, f^^'^ decreases as z ^ 4. In fact, as z becomes larger, also Iuq"] increases and the 
integral in Eq. ( |2.9| ) decreases. However, y^" « —-^0^/2 for z » 1 and since the values of in the E = 3.595 GeV 
NE3 experiment are not large, (j)Q and the corresponding f^^'^ are still non-vanishing at z — > 4. As a consequence, 
the integrals receive a sizeable contribution from the region z > 3 and the NI cross section in the low v region remains 
large. 

Note also that the mixing factor C22 ~ 1 for a; = 0.1 (Fig. 4a) but becomes rather small for a; = 2 (Fig. 4b), sizeably 
reducing the NI part of the cross section in that region. As stated before, without the inclusion of such a factor C22, 
the over-prediction of the theoretical cross section in the l ow ^/ -region would be more pronounced. 

Next, applying the prescription recalled in subsection II C , we make a comparison between GRS and IS cross 
sections for E — 3.595 GeV. Fig. 5a (5b) shows the results for 9 = 30° using the 0o (</'oi) approximation for the 
reduced response, together with the NE3 data. The agreement between the two predictions is good, but not perfect. 
One of the causes is undoubtedly the use of (2.21), which is an approximation for the parametrized, off-shell total 
profile function F, and is of course not intrinsic to the actual ladder summation. Also the agreement with the data is 
good, except for the smallest v. 

One observes that the GRS and IS predictions using only 0o diverge for decreasing v and that the GRS prediction 
is closest to the data. This is shown in Fig. 5a for 9 — 30°, E = 3.595 GeV, but holds in fact for all examined cases. 
Comparison of Figs. (5a) and (5b) moreover shows that FSI for GRS are smaller than for the IS, in particular for 
smaller v. The two observations above can be understood theoretically and have previously been demonstrated 
for simple models. 

Finally, one infers from Fig. 5b that the differences between the order IS and GRS cross sections are noticeably 
reduced when the first order FSI is included in both calculations. A similar comparison is shown in Figs. 5c and 5d 
for 9 = 45° and again the agreement is found to be good after the inclusion of the FSI. 

The separate IS NE and NI parts for E ^ 3.595 GeV and 9 = 30° {9 = 45°) are shown in Fig. 6a (Fig. 6b), where 
we used the same notation adopted in Fig. 3. In this case, the NI part computed with the (j)o approximation for 
the reduced response (thin solid line) stays well below the NE one (thin dashed line) in the low v region, as already 
found in Ref. |3|. Now, in the evaluation of the convolution (2.2), \yo\ becomes rather large as z — > 4 and (/jq rapidly 
decreases (yo — > —\fQ^ as z ^ 1). As a consequence, the IS integrand of Eq. (2.2) is very small in the low v region 
(in contrast to what happens in the GRS case). However, as already discussed in relation with Fig. 5, there the FSI 
contributions in the IS case are sizeable. In fact, now the FSI are given by Eq. ( |2.22D without the term % which 

otherwise would partially cancel the contribution of the large term G^^^ As a result, in the low v region, the IS NI 
cross section calculated at the level of the 0+1 order becomes larger than the NE one, and rather close to the GRS 
NI cross section. 

Another interesting aspect is the influence of the nuclear interaction, chosen to calculate the density matrices. In 
Fig. 7 we display for E = 3.595 GeV, 9 = 30° cross sections computed on the basis of the three previously mentioned 
models of nuclear interaction (AV18UR, AVIS and AV14UR). The calculations have been performed using 0o from 
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Eq. ( 2.15| ) with the density matrices determined directly from the corresponding nuclear wave functions. For identical 
kinematics the predictions for AV18UR and AV14UR can hardly be distinguished, whereas the AV18 cross section is 
slightly different from the other two in the low-v tail: There clearly is only weak dependence on the nuclear interaction. 



IV. CONCLUSIONS 



We have studied inclusive scattering of high-energy electrons from *He, for energy losses below and around the 
quasi-elastic peak and up into the deep inelastic scattering region. The underlying model assumes non-interference 
between nucleonic and sub-nucleonic degrees of freedom, which implies that total nuclear structures function may be 
expressed as a generalized convolution of the structure functions of free nucleons and the one of a nucleus composed 
of point particles. The model is estimated to become gradually imprecise for < 2 — 2.5 GeV^. Structure functions 
for a nucleus of point-particles are computed via the reduced response (f>{q, y) using a relativistic generalization of the 
GRS series, which includes the first order FSI. 

A new element in the development of the latter is an actual calculation of the required single- and two-particle, 
semi-diagonal density matrices, based on accurately computed ''He ground state wave function. The above replaces 
previously used parametrizations of derived density matrices for targets with A > 12. Computed cross sections appear 
to be hardly dependent on the choice of the NN interaction. 

We also exploited a prescription to derive from the GRS series to order 1/q similar terms for the IS. The two 
methods produce cross sections which are rather similar, in particular after the inclusion of the corresponding FSI 
contributions. 

In conclusion our predictions are in good agreement with the NE3 SLAG- Virginia data for all scattering angles, in 
spite of the fact that for 9 < 30° the involved fall below the validity estimate. We also computed cross sections 
for a future JLab experiment at -B = 6 GeV. Its kinematics are largely within the estimated limits of the underlying 
theory and a comparison of theory with data will be more significant for those than is the case for the NE3 data. 
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and for four values of the scattering angles 9 in the (0+l)th order approximation (thick solid lines). Thin dashed (solid) lines 
show the NE (NI) part of the cross sections. Data are from Ref 
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for four values of the scattering angles 9 in the (0+l)th order approximation (thick solid lines). Thin dashed (solid) lines show 
the NE (NI) part of the cross sections. 
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FIG. 3. (a) Contributions to the cross section for E = 3.595 GeV, 6 — 30° from the NE (dashes) and the NI (sohd Unes) 
parts of the nucleon SF. Thin and thick lines are computed from the the 0th and (0+l)th order reduced response ^(\(l\-,yG), 
respectively. The density matrices correspond to the AV18UR nuclear interaction, (b) The same as in (a), but for Q = 45°. 
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factor C22(Q /x , z) (short dashed hne) entering Eq. ( |2.2| ) as function of z for E = 3.595 GeV, 6 — 30° and two cases of x: (a) 
X = 0.1 (corresponding to have = 0.58 GeV^ u = 3.00 GeV) and (b) a; = 2 (Q^ = 2.76 GeV^ u = 0.74 GeV). /(z) is the 
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TABLE I. Selection of cross section ratios He/C, C/Fe. For given angle we give x, and the ratios. Former are from Ref. 
[0, latter from Ref 0. 



6 = 20° 


6 = 25° 


6 = 30° 


X 


Q' 


He/C 


X 




He/C 


X 




He/C 


2.69 


1.44 


0.64 


2.10 


2.07 


0.64 


1.37 


2.51 


0.66 


1.40 


1.34 


0.70 


1.23 


1.88 


0.69 


1.08 


2.35 


1.17 


0.70 


1.18 


0.88 


1.37 


2.51 


0.96 


0.81 


2.12 


1.03 



e = 15° 


6 = 30° 


X 


Q' 


C/Fc 


X 


Q' 


C/Fe 


2.49 


1.05 


0.82 


1.95 


3.38 


0.70 


1.02 


0.97 


1.18 


1.37 


3.09 


0.98 


0.65 


0.91 


0.97 


1.01 


2.79 


1.04 


0.43 


0.83 


1.00 


0.72 


2.43 


1.00 
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TABLE II. Binding energies in MeV of **He calculated with the CHH method using the AV18 and AV18/UIX, and the 
older AV14 and AV14/UVIII, Hamiltonian models. Also listed are the corresponding "exact" GFMC results |33] as well as the 
experimental value. 



Model 




GFMC 



AVIS 

AV18/UIX 
AV14 

AV14/UVIII 
EXP 



24.1(1) 
28.3(1) 
24.2(2) 
28.3(2) 



28.3 



TABLE III. Values of v (in GeV) and (in GeV^) for E = 3.6 GeV ||l|l. The selected values of x correspond approximately 
to the QEP and the lower and upper ends of the curves shown in Fig. 1. 





lower end 


QEP 


upper end 




X 


= 2 


X — 


1 


X = 0.1 


e 


V 


Q' 


V 


Q' 


V Q'^ 


16° 


0.25 


0.93 


0.46 


0.87 


2.00 0.44 


20° 


0.37 


1.40 


0.68 


1.27 


2.50 0.48 


25° 


0.55 


2.06 


0.95 


1.79 


2.80 0.54 


30° 


0.74 


2.76 


1.22 


2.29 


3.00 0.58 


45° 


1.23 


4.99 


1.90 


3.58 


3.31 0.62 



TABLE IV. Values of u (in GeV) and (in GeV^) for E = 6.0 GeV and the values of 9 of the proposed CEBAF experiment. 
The selected values of x correspond approximately to the QEP and the lower and upper ends of the curves shown in Fig. 2. 





lower end 


QEP 


upper end 




X 


= 2 


X = 


1 


X = 0.1 


9 


V 


Q' 


V 




V Q'-' 


23° 


1.21 


4.56 


2.02 


3.80 


5.01 0.94 


30° 


1.80 


6.75 


2.77 


5.20 


5.37 1.01 


45° 


2.90 


10.90 


3.91 


7.34 


5.70 1.07 


60° 


3.69 


13.86 


4.57 


8.58 


5.82 1.09 
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